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FOREWORD 


This  report  was  prepared  by  the  Engineering  Analysis  Division  of 
Computer  Engineering  Associates,  Pasadena,  California,  on  Air  Force 
contract  AF  33 (657)— 8293*  under  Task  Number  821901,  Project  Number  8219. 
The  work  was  administered  by  the  Stability  and  Control  Section,  Control 
Criteria  Branch,  Flight  Control  Division  of  the  Air  Force  Flight  Dynamics 
Laboratory,  Mr.  Robert  L.  Swalm  was  project  engineer  for  the  Laboratory. 

Work  under  the  contract  began  in  March  1962,  and  was  concluded  in 
August  of  1963.  Most  of  the  research  was  carried  out  by  Robert  G. 
Schwendler  and  Jack  H.  Hill,  and  valuable  contributions  were  made  by 
Dr.  Richard  H,  MacNeal,  who  also  directed  the  early  phases  of  the  research 
activity. 

This  is  the  final  report  and  concludes  work  on  contract  AF  33(657)- 
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ABSTRACT 


Tnis  report  contains  an  investigation  of  "mode  interaction"  instabilities 
of  aircraft,  "Mods  interaction"  refers  to  coupling  between  an  elastic  mode 
and  a  rigid  body  mode  of  the  free  system.  Three  ratter  general  airframe 
configurations  are  analyzed  in  detail.  It  is  shown  that  systems  which  tend 
toward  steady  state  divergence  are  particularly  susceptible  to  mode  inter¬ 
action.  Also,  that  aerodynamic  damping  terms  can  have  a  destabilizing 
effect  upon  a  free  system.  The  analyses  show  that  aeroelastlc  systems  which 
possess  no  finite  frequency  elastic  mode  can  be  susceptible  to  a  finite 
frequency  instability.  Computer  studies  of  several  different  airframe 
configurations  are  discussed.  An  appendix  to  the  report  contains  stability 

charts  for  an  aircraft  having  two  rigid  body  modes  and  one  elastic  mode. 
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I.  INTRODUCTION 

The  Investigation  of  refaranea  I,  titled  "Opttmu*  Structural 
Rapresantatlon  tn  Aaroelasttc  Analyses*  was  Initiated  as  a  study  of 
the  affects  of  elastic  sodas  of  an  aaroelasttc  system  on  the  low 
frequency  response  of  the  ays  tan.  The  resulta  of  that  Investigation 
provided  a  straight  forward  nethod  of  representing  a  systas  In  terse 
of  a  few  of  Its  nonsal  coordinates  and  the  "residual  flexibility*  of 
all  higher  sodas  which  proved  to  be  an  accurate  approxlsatlon  of  the 
aeroelasttc  system  for  the  prediction  of  Its  dynamic  behavior  In  the 
frequency  band  from  zero  through  the  frequencies  of  the  normal  modes 
explicitly  Included  In  the  representation* 

This  method  of  structural  representation  was  shown  to  be  valid 
In  all  configurations  studied  except  In  the  case  of  "mode  interaction" 
where  no  conclusions  were  drawn*  "Mode  Interaction"  Is  defined  as  a 
condition  of  potential  or  Incipient  aeroelastlc  Instability  Involving 
one  elastic  mode  and  one  rigid  body  mode  of  the  free  system. 

The  Investigation  reported  here  was  undertaken  to  provide  some 
Insight  Into  the  mechanism  of  "mods  Interaction"  and  to  provide  a 
means  of  predicting  the  susceptibility  of  a  given  configuration  to 
this  phenomenon.  Since  the  "mode  interaction*  phenomenon  Is  definable 
In  terms  of  a  potential  aeroelastlc  Instability*  this  study  was  aimed 
particularly  at  the  prediction  and  understanding  of  the  Instability 
rather  than  at  prediction  of  the  response  at  a  subcrltlcel  spaed* 

Manuscript  released  by  authors  August  196?  for  publication  as  an 
ASD  Technical  Documentary  Report. 

I 


The  prediction  of  critical  velocity  (flutter  speed)  and  the  system 
response  at  a  subcrIHcal  speed  Is  certainly  a  solvable  problem  for  any 
system.  Many  standard  references  on  flutter  analysis  and  aeroelastlct ty 
present  general  methods  which  are  comprehensive  In  their  potential 
application.  Reference  I  presents  equations  of  motion  of  an  aeroelastlc 
system  which  are  an  example  of  a  completely  comprehensive  analysis. 
Equation  (3,69)  of  reference  1  fa 


where  1 


Is  a  diagonal  matrix  of  the  generalized  mass  and  stiffness 
of  the  normal  modes,  -  a,2  £m  1  , 

Is  the  vector  of  normal  coordinate  deflections. 

Is  the  matrix  of  normal  mode  shapes  or  the  "modal  matrix" 
of  physical  coordinates  m.  [*mK]  Is  a  square  matrix, 
[°J  ,S  the  aerodynamic  Influence  coefficient  matrix  In  the 
physical  coordinates  m,  defined  by  ^F^  »  {Vm}  * 

PC]  " ,h*  residual  flexibility  matrix  or  a  matrix  defining 
the  stiffness  properties  of  the  system  which  are  not 
Included  In  the  generalized  stiffnesses  • 

Is  a  vector  of  additional  externally  applied  forces. 


Equation  ( I—  I )  can  provide  a  rigorous  statement  of  the  equations  of 
motion  of  a  system  If  a  1 1  elastic  modes  of  the  system  are  Included  In 
the  K  coordinates  and  at  least  an  excellent  approximation  of  the  system 
equations  when  the  K  coordinates  Include  a  reasonable  nunber  of  the 
lowest  modes  of  the  system. 


2 


Equation  (I-l)  Is  useful  In  providing  an  accurate  prediction  of 
the  dynamic  response  of  a  complicated  system.  Such  predictions,  how¬ 
ever,  do  not  always  give  an  understanding  of  the  basic  mechanism  of 
the  phenomena  being  studied.  Thus,  such  an  equation  can  provide  means 
of  discovering  that  a  given  aircraft  design  has  unsatisfactory  stabil¬ 
ity  and  control  characterlstl cs,  but  unless  the  basic  mechanism  of  the 
Instability  Is  understood  by  the  designer,  he  wl 1 1  be  unable  to  foresee 
which  design  changes  result  In  Improvement.  The  systems  analysed  In 
this  study  were  selected  so  as  to  throw  light  on  the  basic  mechanisms 
of  mode  Interaction.  In  order  to  Illustrate  these  mechanisms  clearly, 

It  was  often  found  useful  to  make  simplifying  assumptions.  Therefore 
the  results  of  this  study  apply  to  more  definite  stability  and  control 
problems  than  does  equation  (I-l),  which  Is  very  general. 

Three  analyses  are  made  In  this  study.  These  analyses  consider 
the  following  systems i 

1.  (Section  II)  Analysis  of  an  airframe  represented  by  two  rigid 
body  modes,  one  elastic  mode,  and  the  "residual  flexibility" 
of  all  hiyher  modes.  The  chief  simplifying  assumption  of  this 
analysis  Is  the  omission  of  all  damping  terms  from  the  equa¬ 
tions  of  motion  of  the  system. 

2.  (Section  III)  Analysis  of  an  airframe  consisting  of  a  simple 
airfoil  flexibly  attached  to  a  rigid  fuselage.  In  addition  to 
the  assumption  that  the  fuselage  Is  rigid.  It  Is  also  assumed 
that  the  airfoil  has  no  mass. 

3.  (Section  3V)  Analysis  of  an  airframe  having  two  rigid  body  modes 
and  one  elastic  mode.  This  analysts  omits  all  consideration  of 
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residual  flexibility  of  the  higher  modes  of  the  system 


Since  this  study  Is  aimed  particularly  at  the  supersonic  and  high 
subsonic  velocity  regime,  aerodynamic  lag  functions  were  omitted  In  all 
analyses. 


k 


II.  ANALYSIS  OF  AN  AIRFRAME  REPRESENTED  BY  TWO  RIGID  BODY  MODES, 

ONE  ELASTIC  MODE,  AND  THE  "RESIDUAL  FLEXIBILITY" 

OF  ALL  HIGHER  MODES 

The  "residual  flexibility"  approximation,  derived  In  reference  I, 

In  generalized  modal  analyses  provides  a  means  of  Including  all  stiffness 
properties  of  a  system  In  the  analyses  while  the  mass  Is  represented  by 
the  generalized  mass  of  a  selected  number  of  lower-frequency  normal  modes. 
Advantages  In  simplicity  of  solution  of  the  conventional  truncated  modal 
approach  are  obtained  when  the  coordinate  velocity  dependent  terms  are 
omitted.  The  omission  of  these  damping  terms  In  the  equations  of  motion 
will  be  detrimental  to  the  accuracy  of  the  predicted  stability  boundaries, 
but  this  analysis  will  be  shown  to  be  useful  by  providing  some  Insight 
Into  the  stability  problems  of  more  complex  structures. 


The  airframe  configuration  considered  In  this  analysis  will  be  a 
general  one  defined  only  by  the  following  parameters, 
m  ■  the  total  mass  of  the  system. 

2 

r  *  the  radius  of  gyration  of  the  system  (then  mr  Is  the 

pitching  mass  moment  of  Inertia  of  the  airframe  about 

Its  c.g.). 
o 

rgm  »  the  generalized  mass  of  the  first  elastic  mode  of  the 
system.  Thus  r^  Is  Just  a  radius  of  gyration  obtained 
by  dividing  the  generalized  mass  of  the  first  elastic 
mode  by  the  total  mass,  and  taking  the  square  root. 


the  undamped  natural  frequency  of  the  first  elastic  mode. 

the  modal  matrix  defining  the  mode  shapes  at  the  h  (plunge) 
coordinates  and  the  6  (pitch)  coordinates  of  the  system 
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for  the  r  modes  Included.  3  modes  are  Included*  the 
zero  frequency  plunge  mode  (mode  l),  the  zero  frequency 
pitch  mode  (mode  2)  and  the  first  elastic  mode  (mode  3). 

the  flexibility  matrix  of  the  complete  free  system 
when  the  zero  frequency  modes  are  restrained  to  zero 
displacement,  partitioned  by  h  and  0  coordinates. 

the  "residual  flexibility"  matrix  of  all  elastic  modes 
higher  than  the  first,  partitioned  by  h  and  0  coordinates. 
This  matrix  actually  Is  not  needed  to  define  the  system 
since  (from  a  special  case  of  equation  (3.21)  of  refer¬ 
ence  l)  It  Is  determined  from  system  parameters  listed 
above  by  the  equation 


>- 
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v 

kl 

r 1  i 
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X!  1 
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where  Xjj  Is  the  flexibility  matrix  with  the  zero  fre¬ 
quency  modes  restrained. 


The  aerodynamic  forces  Included  In  this  analysis  will  be  defined  by 
the  following  partitioned  matrix  equation 


M 

1 

I 

O 

o 

:x 

o. 

< 

h 

KJ 

-  -  q 

0  1  0 

1 

— 

0 

where  q  Is  the  dynamic  pressure  ^  P  /  , 


The  equations  of  motion  of  the  aeroelastlc  system  can  be  written  by 

2 

a  process  of  partitioning  equation  (i-l).  The  quantity  -  w  In  equation 

p 

(I-l)  Is  replaced  by  s  In  equation  (II-2),  since  we  look  for  solutions  of 
the  equations  of  motion  having  time  dependence  of  the  form  es^,  where  s 
may  be  complex. 
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In  order  to  shorten  subsequent  matrix  equations  we  shall  define  an 
"aeroelastle  aerodynamic  Influence  coefficient  matrix"  ^hfi]  as  *°H0WS 

■  {[J]  -  i  [«»]  [*£]}"  [’he]  (“-« 

It  Is  well  to  note  that  Is  a  function  of  the  dynamic  pressure 

(q)  and  the  number  of  elastic  modes  explicitly  Included  In  the  analysis 
(In  this  case  I  mode). 

The  equations  of  motion  of  equation  ( II— 3)  may  he  further  simplified 
by  recognizing  that  the  modal  colimn  {*.■}  ,  the  pitching  slope  deflec¬ 
tions  In  the  rigid  plunge  mode.  Is  equal  to  zero,  and  that  the  modal  col¬ 
umns  M  and  M  ,  the  plunge  deflections  In  the  rigid  plunge  mode 
and  the  pitch  slopes  In  the  rigid  pitch  mode,  are  equal  to 
column  of  ones,  due  to  normal Izatlon.  The  equations  of  motion  now  may 
be  written  In  the  form 

-  c-}T  <  pj  {♦.’}  1  [V 

0.  0  »rV  -  {V}T  ,  [?„,]{)}  -{»h8}T’[«h.]{».J}  •  • 

'  {*h5}  ’  [?h«]  {'}  -r.S(52,“.2)-{*h5}  ’  |?m]  {‘4 

(13-5) 

In  an  effort  to  obtain  a  solution  In  terms  of  parameters  for  which 
engineers  are  more  likely  to  have  Intuitive  judgment  than  for  the  general¬ 
ized  modal  parameters  of  equation  (II-5),  we  wl 1 1  Introduce  the  concepts 
of  steady-state  elastic  divergence  and  short  period  mode  frequency* 

Conventionally  the  short  period  mode  frequency  Is  the  pitching  fre¬ 
quency  exhibited  by  the  rigid  system  when  aerodynamic  forces  are  Included. 


0 


For  this  analysis  we  will  define  the  "aeroelastlc  short  period  mode 
frequency"  («Q*)  to  be  similar  to  the  conventional  definition  except 
that  the  system.  Instead  of  being  rigid,  will  Include  the  "residual 
flexibility"  of  all  elastic  modes  higher  than  the  first.  This  fre¬ 
quency  can  be  obtained  from  equation  (II-5)  by  letting  be  zero, 
then  assuming  harmonic  motion. 


( 1 1- 6) 


Steady-state  elastic  divergence  Is  conventionally  defined  for  a 
supported  system  as  that  aeroelastlc  Instability  which  can  occur  at 
zero  frequency  and  Is  therefore  Independent  of  the  mass  of  the  system. 
This  simple  definition  cannot  be  applied  to  a  free  system  because  In 
this  case! 

1.  Other  potential  zero-f requency  Instabilities  exist  which 
are  not  aeroelastlc  In  nature. 

2,  The  mass  distribution  of  a  free  system  must  be  considered. 

The  only  reasonable  definition  for  steady-state  elastic  divergence  of 
a  free  system  Is  available  through  the  description  of  the  elasticity 
of  the  str  -stupe  when  restraints  are  placed  on  the  zero-f  requency 
normal  modes  of  the  system.  In  this  case  steady-state  divergence  can 


{H  -  Ml  M} 


be  defined  by  the  singularity  of  the  matrix 
where  q  |q  J  is  the  matrix  of  aerodynamic  Influence  coefficients  and 
Is  the  deflection  Influence  coefficient  matrix  of  the  structure 
when  the  zero-frequency  modes  are  restrained.  In  this  analysis  elastic 
divergence  Is  determined  by  the  vanishing  of  the  determinant 


H  - 1  MM 
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which  may  be  expanded  as  follows  In  the  parameter  q 


>]  ■  ’  Ml>] 


■  '  "  0,q  +  D2q  - 


where:  D|  Is  the  sum  of  the  first  symmetric  minors  of  j^hej[*6hj 
(sum  of  the  diagonal  elements), 

Dg  is  the  sum  of  the  second  symmetric  minors  of 
0^  Is  the  sum  of  the  h*h  symmetric  minor  of  £xgh  . 

Then  defining  the  "aeroelastlc  Index"  D  by 


( 1 1-7) 


D.q  -  D0q  + 


(II-8) 


from  equation  ( II-7 ) J  D  wl 1 1  equal  0  when  q  »  0  and  D  wi 1 1  equal  I  when 
q  Is  the  dynamic  pressure  at  divergence. 


Steady-state  divergence  of  the  free  system  can  also  be  determined 
from  equation  ( II-5)  by  letting  ?!  ■  ^  «  s  •  0.  Then 

mre2  -{*h5}  qD  phe]  {V}  -  0  (n“9) 

where  qQ  Is  the  dynamic  pressure  at  divergence.  We  now  note  that  the 
quantity 

e  e 

exhibits  the  same  known  dependence  on  q  as  Dj  that  Is  both  functions: 

1.  Equal  I  when  q  *=  0. 

2.  Equal  0  when  q  *  the  dynamic  pressure  at  divergence, 

3.  Are  equal  to  q  QX^g  when  applied  to  a  system  containing 
only  I  elastic  degree  of  freedom,  1  pair  of  aerodynamic 
coordl nates. 

It  is  then  postulated  that 

5  ■  {V? q  W  {n5}  '“-'d 

e  e 


10 


A  rigorous  proof  of  equation  (ll-ll)  was  not  obtained  and  Its  validity 
will  be  tried  by  example  problems. 


The  characteristic  equation  obtained  from  the  determinant  of  the 
matrix  coefficients  of  equation  ( II— 5)  Is 


s2  [s'*  +  As2  ♦  b]  »  0 


(II- 1 2) 


where 


■  {v>  MT  ■>  W  {■}} 


For  a  configuration  where  the  slope  deflections  "sensed"  by  the 
aerodynamic  representat Ion  are  all  equal  In  the  first  elastic  mode,* 
then  the  last  term  In  equation  ( II— 1 2)  for  the  factor  B  Is  zero,  since 
through  normalization  can  be  made  equal  to  <£l  .  For  the 

configuration  considered  In  Section  IV  (since  only  1  aerodynamic  cell 
Is  Included)  this  term  Is  zero  and  probably  for  a  rather  large  range 
of  configurations  this  term  has  little  Importance.  Therefore  because 
this  term  Is  of  doubtful  significance  and  because  Its  omission  greatly 
simplifies  the  results  of 'this  analysis  It  will  be  assumed  that 


(II- 1 5) 


*An  example  of  this  configuration  Is  one  containing  a  rigid  lifting 
surface. 


1 1 


Then  by  substitution  of  equattons  (II-6),  ( II-- It)  and  (11-13)  Into 
equation  ( II— 12),  the  characterlstl c  equation  con  be  written 


s2  |>  +  As2  ♦  b]  «  0, 

where 

2  *2 
A  -  0>e  (  I  -  B)  +  tt 

„  2 
B  »  to  to 
o  e 


( xi-i  U) 


—  2 

In  equation  (II-I1+)  wl th  ( 1  -  B)  2  0,  a  wl  1 1  be  real  and  negative 

p  p 

unless  A'  -  ItB  Is  negative*  If  Ac  -  I4B  Is  negative,  then  two  of  the 
roots  a  of  equation  (II-I1+)  will  have  a  positive  real  part.  The 
system  Is  therefore  on  the  verge  of  Instability  when 

A2  -  1|B. 

Then  marginal  stability  Is  given  by 

* 

(,  ,  cii-15) 

« * 

and  since  both  D  and  are  always  continuous  functions  of  the  dynamic 

6 

pressure  beginning  at  the  origin,  the  lowest  dynamic  pressure  which  will 

satisfy  equation  (II-I5)  Is  given  by 

* 

-  1  -  Vb7  (11-16) 

0 

Both  sides  of  equation  (II-I6)  are  functions  of  the  dynamic  pressure. 

The  dynamic  pressure  at  flutter  Is  defined  as  the  Intersection  of  these 
2  functions. 


The  first  observation  we  may  make  from  the  results  of  this  analysis 


Is  that  a  strong  relationship  does  exist  between  mode  Interaction,  or 


flutter  resulting  from  coupling  of  the  first  elastic  mode  and  the  rigid 
body  modes,  and  elastic  divergence.  This  analysts  shows  that  for  free 
systems  which  fit  the  assumptions  made,  classical  "elastic  divergence” 
will  not  occur  because  a  flutter  Instability  will  always  exist  at  a 
velocity  lower  than  that  correspond! ng  to  divergence. 

The  second  observation  we  may  make  from  equation  (II-I6)  Is  that 
when  0,  the  "aeroelastlc  Index",  Is  negative  the  frequency  ratio  Is 
complex  and  therefore  flutter  will  not  occur.  The  significance  of  a 
negative  D  Is  that  the  system  Is  losing  aerodynamic  effectiveness  (at 
a  given  rigid  body  mode  pitch  angle,  the  aeroelastlc  deflections  of  the 
structure  reduce  the  total  lift  force).  This  observation  provides  the 
following  useful  qualitative  criteria  for  the  susceptibility  of  a  given 
configuration  to  mode  Interaction! 

1.  If  a  system  loses  aerodynamic  effectiveness  as  velocity  Is 
Increased,  then  mode  Interaction  will  not  occur, 

2.  If  a  system  Increases  In  aerodynamic  effectiveness  as  velocity 
Is  Increased,  then  mode  Interaction  will  probably  occur;  the 
system  will  be  unstable  at  a  velocity  less  than  that  predicted 
for  steady-state  divergence.  At  the  velocity  of  Instability 
the  ratio  of  the  uncoupled  short  period  mode  frequency  to  the 
first  elastic  mode  frequency  will  be  less  than  I. 

COMPARISON  WITH  KNOW  SOLUTIONS 

A  number  of  aircraft  were  studied  on  the  CEA  analog  computer  In 
the  course  of  this  project  and  that  reported  In  reference  I.  Most 
configurations  studied  exhibited  a  flutter  Instability  Involving  the 
rigid  body  modes  and  the  first  elastic  mode,  and  In  all  but  one  of 
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these  cases  the  system  tended  to  diverge  at  some  higher  speed.  One  con¬ 
figuration,  (configuration  3  of  reference  1)  a  swept  wing  airplane,  lost 
aerodynamic  effectiveness  as  velocity  Increased.  This  configuration 
showed  no  tendency  to  flutter  even  though  the  short  period  mode  frequency 
virtually  coincided  with  the  elastic  mode  frequency.  The  one  configura¬ 
tion  which  did  not  tend  to  diverge  but  dtd  exhibit  mode  Interaction  was 
one  where  finite  aerodynamic  damping  wl  1 1  be  shown  to  be  a  necessary 
parameter  for  flutter  to  occur  and  therefore  could  not  be  predicted  by 
this  analysis  which  Ignored  all  damping  terms. 

These  results  In  general  coincide  with  the  qualitative  criteria 
drawn  from  this  analysts}  however.  In  most  cases  the  aerodynamic  damp¬ 
ing  appeared  to  have  a  reasonably  large  effect  on  the  quantitative  value 
of  the  flutter  speed,  precluding  accurate  prediction  by  this  analysis. 

Configuration  4  of  reference  I,  the  delta  winged  airplane,  showed 
very  rapid  variations  In  damping  at  flutter.  This  phenomenon  usually 
Indicates  that  aerodynamic  damping  terms  have  little  Influence  on  the 
Instability  because  they  do  not  vary  rapidly  with  speed  or  A.C.  location. 
The  criterion  of  equation  (II-I6)  was  applied  to  configuration  4  because 
damping  did  not  appear  to  be  a  controlling  parameter  In  this  case.  This 
criterion  was  the  only  one  produced  In  this  study  which  was  In  a  form 
applicable  to  this  "plate-llke"  structure. 

DESCRIPTION  OF  CONFIGURATION  4 

The  aircraft  used  in  this  comparison  Is  Identical  to  configuration 
4  of  reference  I.  The  geometry,  mass  distribution  and  structural  param¬ 
eters  are  repeated  in  this  report  In  Figures  I  through  5»  The  aero- 


dynamic  forces  were  represented  as  follows! 

1.  For  the  purposes  of  describing  aerodynamic  forces,  the  wing 
was  divided  Into  three  strips  shown  In  Figure  6.  The  strips 
are  assured  to  be  rigid  planes  whose  deflection  Is  defined  by 
the  plunge  deflections  at  the  \/k  and  3A  chord  coordinates 
shown  In  Figure  6.  The  aerodynamic  lift  and  moment  on  each 
strip  are  rigidly  beamed  to  these  same  i/k  and  3/U  chord  co¬ 
ordinates  of  the  elastic  structure. 

2.  The  aerodynamic  center  of  each  strip  is  located  on  the  mean 
chord  of  the  strip  and  aft  of  the  effective  leading  edge  of 
the  strip  a  distance  xc  (c  ■  mean  chord  length  of  strip). 

x  was  varied  In  the  study. 

3.  The  aerodynamic  lift  force  on  a  strip  Is  given  by 


where 

6  *  |  ]  Is  the  pitching  slope  of  the  strip 

(positive  nose  up), 

z  Is  the  plunge  deflection  at  the  A.C.  (positive  up), 
and  Zj  are  the  plunge  deflections  at  the  l/U 
chord  and  3/k  chord  points  of  the  elastic  structure 
(positive  up). 

Thus 

£ 

z-Zi(|-2x)  -z?(^-2x). 

*This  equation  Is  given  Incorrectly  In  reference  I. 
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4»  The  aerodynamic  moment  about  the  A.C.  Is  given  by 

2 

Ik  •  «  ^  P  V2  S  ^  (positive  nose  up). 

5.  The  basic  flight  condition  used  tn  the  following  numerical 
study  Is  described  by  the  parameters* 

Velocity  ■  1655  mph 

Altitude  »  40,000  feet 

/  2 

Dynamic  Pressure  »  qQ  ■  11.9  lb. /In. 

Lift  Curve  Slope  ■  *  5»0  per  rad. 

a 

The  dynamic  pressure  Is  varied  In  the  numerical  study  and  Is 
expressed  as  fraction  of  the  basic  value  given  above  (qQ). 

NUMERICAL  COMPARISON 

This  configuration  was  simulated  on  the  CEA  passive  analog  computer 
In  the  project  reported  In  reference  I.  At  that  time  It  was  observed  to 
exhibit  the  following  unusual  aeroelastlc  properties* 

1.  Three  distinct  flutter  Instabilities  were  observed  for  various 
combinations  of  dynamic  pressure  and  A.C.  location.  Two  of 
these  Instabilities  were  Identified  as  conventional  flutter 
phenomena  Involving  the  coupling  of  2  elastic  modes.  The  third 
Instability  was  the  result  of  coupling  of  elastic  modes  and 
rigid  body  modes.  It  Is  the  latter  case  which  wll  I  be  consider¬ 
ed  In  this  comparison. 

2.  The  system  was  observed  to  have  a  much  smaller  stable  region 
when  the  system  was  represented  by  one  elastic  mode  than  when 
the  higher  modes  were  Included  In  the  representation. 
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3«  The  stability  of  the  system  was  sensitive  to  small  variations 


In  flexibility. 

The  application  of  equation  (II-I6)  to  configuration  Involved 
two  major  steps t 


I.  The  determination  of 


for  various  A.C.  locations  and 


values  of  -3-  .  This  calculation  was  accomplished  using 
% 

equations  (II-6),  (II-U)  and  the  modal  properties  of  the 
system  presented  In  reference  I. 

2.  The  calculation  of  D  from  equation  (II-8)  and  the  flexi¬ 
bility  matrix  of  the  system  given  In  reference  1. 

The  results  of  these  calculations  are  presented  graphically  In  Figures 
7  through  II  for  values  of  x  of  .25.  .30*  .35*  .375  and  .I4D.  Also 
Included  on  Figures  7  through  II  are  curves  of  the  dimensionless  short 


period  mode  frequency  of  the  rigid  system 


to 


and  the  curve  of  I 


-V? 


when  the  elasticity  of  the  system  Is  represented  by  only  the  generaliz¬ 
ed  flexibility  of  the  first  elastic  mode.  The  Intersection  of  the  latter 
two  curves  defines  the  dynamic  pressure  ratio  °t  marginal  stability 

when  the  "residual  flexibility"  of  all  elastic  modes  higher  than  the 
first  Is  Ignored.  The  Intersection  of  the  upper  curves  defines  the 
stability  boundary  for  the  system  when  all  flexibility  of  the  system 
Is  Included. 


The  stability  boundaries  calculated  from  equation  (II-I6)  are  pre¬ 
sented  in  Figure  12  along  with  the  corresponding  boundary  observed  In 
reference  I. 
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Curve  A  of  Figure  12  Is  the  low  frequency  stability  boundary  of 
reference  I  (dynamic  pressure  ratio  has  been  substituted  for  the 
equivalent  parameter  "flexibility  factor"  of  reference  I  for  convenience). 
Curve  B  of  Figure  12  Is  the  stability  boundary  calculated  from  equation 
( II— 1 6)  Including  the  residual  flexibility  of  all  modes  of  the  system. 
Curves  A  and  B  are  In  substantial  agreement.  The  discrepancy  between 
curves  A  and  B  can  be  attributed  to  any  one  or  all  of  the  following 
reasons. 

1.  The  approximations  made  In  this  analysis. 

2.  The  omission  of  all  damping  terms  from  this  analysis. 

3.  Experimental  error  In  the  determination  of  system  flexibility. 
Olscrepanctes  were  also  observed  on  the  analog  computer  between 
the  modal  simulation  and  the  "exact"  representation.  The  sta¬ 
bility  boundaries  were  recorded  for  the  "exact"  representation. 

Curve  C  of  Figure  12  Is  the  stability  boundary,  determined  from 
equation  (II-|6),  when  the  flexibility  of  the  system  Is  represented  by 
the  generalized  flexibility  of  only  the  first  elastic  mode.  The  stable 
region  above  curve  B  Is  much  larger  than  that  above  curve  C.  This  fact 
agrees  with  observations  of  the  analog  computer  analysis  of  reference  I. 

It  is  not  claimed  that  this  comparison  proves  the  validity  of 
equation  (II-I6).  However  It  Is  known  that  equation  (II-I6)  Is  valid 
for  very  simple  systems  and  this  comparison  does  show  that  It  Is  useful 
for  qualitative  evaluation  of  a  complicated  system.  The  conclusion  Is 
made  that  the  criteria  of  equation  (II-I6)  Is  useful  In  the  prediction 
of  the  susceptibility  of  a  configuration  to  the  "mode  Interaction" 
phenomenon. 
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Notet  Numbered  Values  Show  Spar 
Sending  Moment  of  Inertia 

E  -  I07  lb./tn.2 

l/l  Varies  Linearly  Between 
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Figure  2,  Spar  Stiffness  of  Configuration  1+ 
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Figure  lu  Torque  Boxes  of  Configuration  4 
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Dimensionless  Dynamic  Pressure 
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Figure  10. 
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Figure  12 
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III.  ANALYSIS  OF  AN  AIRFRAME  CONSISTING  OF  A  SIMPLE 
AIRFOIL  FLEXIBLY  ATTACHED  TO  A  RIGID  FUSELAGE 

In  order  to  place  frequency  coalescence  methods  In  proper  perspective. 
It  Is  useful  to  analyze  a  basic  aeroelastlc  system  consisting  of  a  simple 
rigid  lifting  surface  flexibly  mounted  to  a  rigid  fuselage  (see  Figure  I3)« 
Since  we  are  not  concerned  here  with  conventional  binary  flutter,  we  shall 
assume  that  the  inertia  of  the  lifting  surface  Is  entirely  negligible. 

Of  course,  by  Ignoring  the  Inertia  of  the  lifting  surface,  we  leave  out 
of  account  Instabilities  which  occur  because  of  coupling  between  the  short 
period  mode  and  a  low  frequency  elastic  mode.  This  type  of  Instability 
will  be  discussed  In  Section  IV  of  the  report. 

Clearly,  If  the  system  possesses  Instabilities,  they  are  not  of  the 
type  which  can  be  predicted  by  frequency  coalescence  methods,  since  there 
Is  no  finite  frequency  elastic  mode  which  can  coalesce  with  the  short 
period  mode.  We  will  demonstrate  that  this  system  can  exhibit  three 
distinct  types  of  Instability, 

The  oonf I guratlon  of  the  system  can  be  specified  by  four  coordinates* 

z  ■  vertical  deflection  of  the  center  of  mass,  positive  up, 

8  ■  pitch  angle  of  the  fuselege,  positive  nose  up, 

z  -  vertical  deflection  of  the  aerodynamic  center,  positive  up, 

8 

8  “  pitch  angle  of  the  lifting  surface,  positive  nose  up. 

d 

All  coordinates  are  measured  relative  to  an  Inertial  reference  system. 

The  equations  of  motion  are 

mz  -  K(za  -  z  -  xe8  -  x8e  8a)  -  L 

I®  xe  K<2a  ”  2  "  xe®  "  *ae  ®a) 
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o  -  -  K(za  -  Z  -  xe0  -  xa6  ea)  +  L 

0  *  xae  K(za  "  z  “  xe6  "  xae  ®a}  '  k(°a  "  4  M 

The  lift  and  moment  at  the  aerodynamic  center  are 


i  c  e 

L  .  a  s  Ic  e  -  c  ~+c  _ - 

H  >  L  a  L  V  L  V 
a  a  q 


M  “  q  S  -  C 


c2$. 


q  ■*  •£  P  V2 

(the  geometrical  quantities  xa,  xe,  xflg  are  defined  In  Figure  13). 
Introducing  matrix  notation,  these  equations  may  be  written* 
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Introduce  c  as  a  unit  of  length,  ^  as  a  unit  of  time,  and  PSc  as 
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a  unit  of  mass.  Then  the  following  quantities  are  non-dimensional* 
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In  terms  of  these  quantities  the  equations  of  motion  become 
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and  the  characteristic  equation  Is  obtained  by  taking  the  determinant 
of  the  matrix  of  coefficients.  Laplace's  method  of  expanding  a  deter¬ 
minant  (reference  5)  Is  convenient  and  yields  the  following  characteristic 
polynomial t 

S2[a^+  a?53+  OgS2  +  a|5+  aj  -  0, 


where 
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Before  deriving  any  conclusions  concerning  the  stability  of  the 
system,  It  Is  useful  to  Introduce  a  fictitious  frequency,  toQ  ,  which 
would  be  the  undamped  short  period  frequency  If  the  airframe  were  rigid. 
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The  usefulness  of  this  parameter  depends  upon  the  fact  that  given  ci)Q,  the 
velocity  V  Is  determined  and  vice  versa,  and  upon  the  fact  that  stabil¬ 
ity  criteria  are  more  conveniently  expressed  In  terms  of  frequency  ratios 
than  In  terms  of  velocity.  The  non-dimensional  frequency  corresponding 
to  6)q  Is  simply 


In  the  characteristic  equation  replace  by 
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and  divide  the  equation  by  r^  os  In  this  way  we  obtain 
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Obviously,  the  characteristic  polynomial  has  a  pair  of  zero  roots. 
These  roots  merely  express  the  fact  that  all  altitudes  and  all  directions 
of  flight  are  equivalent.  There  are  several  Interesting  special  cases. 


A,  Suppose  the  torsional  stiffness  Is  Infinite  (5>e  »  °°). 
In  this  case  the  characteristic  polynomial  Is 

s^  +  b„  s2  +  b  S  +  b  «  0, 
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For  a  given  airframe.  It  Is  convenient  to  take  the  frequency  ratio 


“o 

—  as  a  measure  of  velocity.  Since  w.  Is  known  and  Is  Independent  of 
CO  z 

2  “o 

velocity,  —  Is  proportional  to  velocity  V.  The  value  of  V  at  which 
U2 

the  airframe  Is  neutrally  stable  Is  of  particular  Interest,  and  can  be 
determined  as  follows.  Suppose  the  airframe  Is  neutrally  stable,  then 
there  Is  at  least  one  root  of  the  form  s  *  IE  where  to  Is  real.  It 
follows  Immediately  from  the  characteristic  polynomial  that  if  u  /  Os 
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This  condition  determines  the  value  of 
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In  many  cases,  especially  In  cases  where  the  surface  Is  mounted 
well  aft  of  the  center  of  mass  of  the  airframe,  the  effects  of  and 


C  (pitch  damping  coefficients)  will  be  negligible  compared  to  the 

mq  C.  C_ 

L  m 

effects  of  C.  .  In  such  a  case  we  can  Ignore  the  ratios  ——1  ,  , 

L  CL  L 
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and  an  especially  simple  expression  for  the  frequency  ratio 
neutral  stability)  Is  obtained* 
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This  criterion  can  be  given  a  simple  physical  characterization.  Imagine 
the  aerodynamic  center  to  be  rigidly  constrained,  the  airframe  being 
otherwise  free  of  forces  and  constraints,  then  the  natural  frequency 
of  the  system  Is 


(III- 1 6) 


Therefore,  when  the  short  period  frequency  coQ  (computed  as  though 
the  airframe  were  rigid)  equals  the  frequency  to*,  neutral  stability 
occurs.  The  frequency  of  the  neutrally  stable  oscillations  of  the  alr- 
2  ^o  2 

frame  Is  simply  u  ■  ,  the  short  period  frequency  computed  by 

assumtng  the  airframe  to  be  rigid.  Although  the  presence  of  pitch  damp¬ 
ing  due  to  CL  and  will  modify  these  conclusions  slightly,  the  basic 

q  q 

mechanism  of  the  Instability  for  practical  cases  Is  revealed  by  (III-I5) 
and  (III- 1 6). 

Notice  that  the  frequency  ratio  at  neutral  stability  can  be  express¬ 
ed  In  terms  of  four  dimensionless  parameters* 
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and  If  damping  terms  are  small. 


only  a  single  parameter  Is  Involved, 


r2 

ac 


» 


D.  Infinite  plunge  stiffness  («z  ■  oo  ) . 

The  characteristic  polynomial  Is  again  of  third  order* 
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Suppose  s  ■  IS,  U  real.  Then 

-  c,  tS^  +  c.  <S  -  0, 
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Obviously,  to  ■  0  Implies  c  *  0,  and  this  Is  not  possible  unless  —  -  0. 

o  ®0 

This  last  condition  Implies  that  the  surface  is  connected  to  the  airframe 

by  means  of  a  torsional  spring  of  zero  stiffness  or  else  the  velocity  Is 

Inf  Ini  te,  I  f  to  /  0  then 
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N°w  xac  ts  negative  It  the  airframe  has  a  positive  static  margin.  There¬ 

of 

fore  If  the  elastic  axis  Is  aft  of  the  aerodynamic  center,  - Is 

xac 

positive,  and  Instability  Is  certainly  possible.  It  appears  that  the 
existence  of  an  Instability  Is  closely  connected  with  the  divergence 
properties  of  the  lifting  surface.  Thus,  If  -  "x  Is  positive,  the 
surface  Is  capable  of  divergence  when  the  fuselage  of  the  airframe  Is 
restrained.  The  Instability  predicted  by  (III-22)  Is  as  close  as  the 
free  system  can  come  to  divergence.  The  Instability  Is  not  ordinary 
divergence  since  the  frequency  of  the  Instability  Is  greater  than  zero. 


It  Is  not  generally  possible  to  Ignore  the  pitch  damping  terms  C  , 

q 

In  equation  (XXI-22),  since  the  character  of  the  Instability  can  be 

q 

drostical ly  modi f I ed  by  these  terms.  In  fact,  there  Is  one  Interesting 
case  where  the  pitch  damping  terms  can  produce  instability  at  a  very 
low  airspeed.  Suppose  we  position  the  aerodynamic  center  with  respect 
to  the  elastic  axis  In  such  a  way  that 


(I 11-23) 
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We  con  then  position  the  lifting  surface  with  respect  to  the  center  of 
gravity  of  the  airframe  so  that 

Cl  Cm 

<ra-2« 

a  Ua 

Then  it  Is  obvious  from  the  expressions  given  for  the  coefficients  of 
the  characteristic  polynomial  that 
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and  therefore  the  system  can  be  neutrally  stable  for  all  values  of 


(e.g.  the  system  Is  neutrally  stable  at  every  airspeed  from  0  to°o). 
Notice  that  Is  crucial  to  this  phenomenon.  If  ■  0  It  Is  not 

possible  to  satisfy  either  of  equations  (III-23)  or  (III-21*).  It  would, 
of  course,  be  rare  to  encounter  an  airframe  which  satisfied  the  conditions 
of  equations  ( III— 23) •  (111-21*)  exactly.  However  the  example  does  show 
that  the  effects  of  pitch  damping  can  be  very  Important,  Note  also  that 
structural  damping  (dashpot  In  parallel  with  the  torsional  spring  of 
Figure  I j)  would  tend  to  eliminate  this  Instability  at  low  airspeeds. 


The  frequency  of  the  Instability  Is  apparent  from  equation  ( III— 25) • 


In  the  general  case,  where  both  Ug  and  &>2  are  finite,  the  substitu¬ 
tion  s  -  115  Into  the  characteristic  equation  ( III— 2*)  yields 
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A,  -  A„  S2  +  A  -  0 
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where  the  A's  are  given  by  equations  (III-5),  (III-6),  (III-7),  (III-8), 
(III-9).  If  S  /  0,  these  two  equations  can  be  combined  Into 
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Happily,  If  this  equation  Is  rearranged.  It  becomes  linear  In 
and  we  ftnd 
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The  frequency  of  oscillation  Is  obtained  from  Q  ■  -j—  : 
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The  frequency  ratio  —  at  neutral  stability  Is  determined  by  the 


following  set  of  seven  non-dimensional  parameters: 
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Our  conclusions  regarding  airframes  which  consist  of  a  single  light 
lifting  surface  elastically  mounted  to  a  rigid  fuselage  may  be  summarized 
as  fol lows i 

A,  The  system  may  be  dynamically  unstable  even  though  It  Is 
statl ca 1 1 y  stable. 

B.  The  velocity  at  which  the  airframe  becomes  neutrally  stable 

“e 

can  be  determined  from  the  frequency  ratio  ~  ,  which  can 


be  expressed  In  terms  of  the  seven  non-dimensional  parameters 

C 

_  2  mq 

listed  above,  The  parameter  UQ  Is  usually  quite  small 


a 

and  Is  less  Important  than  the  other  parameters. 

C.  If  the  torsional  stiffness  between  the  lifting  surface  and  the 

fuselage  is  infinite,  an  instability  will  occur  when  the  short 

♦ 

period  frequency  «0  equals  the  frequency  o>  .  The  latter  fre¬ 
quency  is  simply  the  natural  frequency  of  the  fuselage  on  the 


12 


plunge  stiffness  when  the  displacement  of  the  lifting  surface 
at  the  aerodynamic  center  Is  restrained.  This  simple  criterion 
Is  strictly  correct  only  If  the  effects  of  aerodynamic  pitch 
damping  are  Ignored!  however,  the  criterion  Is  not  sensitive 
to  pitch  damping  for  practical  airframes. 

0.  If  the  plunge  stiffness  between  the  lifting  surface  and  the 
fuselage  Is  Infinite,  el ther  of  two  types  of  Instabl 1 1 ty  may 
occur.  If  the  aerodynamic  center  Is  forward  of  the  elastic 
axis  (x  <  0),  an  Instability  occurs  which  Is  connected  with 
the  possibility  of  lifting  surface  divergence.  In  addition, 
by  careful  positioning  of  the  elastic  axis  with  respect  to  the 
aerodynamic  center  and  the  center  of  gravity.  It  Is  possible 
to  obtain  an  airframe  which  Is  neutrally  stable  at  all  air¬ 
speeds.  This  type  of  Instability  can  only  be  produced  If  CL 
Is  different  from  zero. 

E,  None  of  the  above  Instabilities  can  be  predicted  on  the  basis 
of  frequency  coalescence  methods.  We  can  describe  the  Insta¬ 
bilities  as  being  due  to  aerodynamic  coupling  between  zero 
frequency  modes  (rigid  motions)  and  Infinite-frequency  modes, 
since  the  structure  In  vacuo  possesses  only  these  types  of 
modes. 

F.  The  conclusions  of  this  section  should  apply  to  airframes  with 
a  simple,  light,  lifting  surface  and  a  fairly  rigid  fuselage. 

A  more  realistic  type  of  airframe  Is  analyzed  In  Section  IV. 
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5  gure  13.  Airframe  Consisting  of  a  Simple  Airfoil  Flexibly 
Attached  to  a  Rigid  Fuselage 


IV.  ANALYSIS  OF  AN  AIRFRAME  HAVING  TWO  RIGID 
BODY  MODES  AND  ONE  ELASTIC  MODE 


The  configuration  shown  In  Figure  |j+  Is  In  some  ways  more  general 
than  the  configuration  analyzed  In  Section  III  and  Is  also  a  configura¬ 
tion  which  frequently  occurs  In  practtce.  A  detailed  analysis  should 
therefore  be  useful.  It  Is  convenient  to  give  the  equations  of  motion 
In  modal  form.  There  are  two  rigid  body  modes  and  one  elastic  mode. 

If  the  rigid  mode  shapes  are  properly  normalized,  the  equations  of 
motion  of  such  a  system  can  always  be  written  In  the  form 
mz  *  / P(x,  y)<3S 
Id’  *  /xp(x.  y)dS 

n>e(€  +  «e2  £)  “  A(x,  y)p(x,  y)dS 

where 

m  »  the  total  mass  of  the  system 
I  ■  the  total  pitching  Inertia  about  the  center  of  mass 
mg  ■  the  generalized  mass  of  the  elastic  mode 
p(x,  y)  -  net  vertical  pressure  on  the  element  of  area 
dS  at  the  point  (x,  y) 

*(x,  y)  ■  mode  shape  of  the  elastic  mode 

If  £  ■  0,  then  the  two  normal  coordinates  r  and  8  have  a  simple  geomet¬ 
rical  significance,  z  being  the  vertical  displacement  of  the  center  of 
mass  and  8  the  pitch  angle  of  the  fuselage.  The  geometrical  significance 
of  £  depends  upon  the  way  In  which  the  mode  shape  #(x,  y)  Is  normalized. 
We  shall  assume  that  a  single  rigid  lifting  surface  Is  mounted  at  some 
point  on  the  fuselage  as  shown  In  Figure  li+.  Then  .the  vertical  displace¬ 
ment  of  any  point  of  the  surface  Is 

z  +  x8  +  4>(x,  y)£ 

1*5 


and  the  slope  of  the  lifting  surface  Is 


©  + 


d$ 

3x 


obtained  by  taking  the  partial  derivative  of  the  vertical  deflection 

d$ 

with  respect  to  x.  Since  the  surface  Is  assumed  to  be  rtgld,  -5-  must 

ox 

be  a  constant  over  the  surface.  If  we  normalize  the  elastic  mode  shape 
3$ 

by  setting  ^  ■  I,  then  §  has  a  simple  geometrical  significance.  When 

6  -  0,  4  Is  simply  the  pitch  angle  of  the  lifting  surface.  It  follows 
3$ 

from  •£-  m  |  that  4>(x,  y)  -  x  +  constant.  It  is  convenient  to  express 
this  constant  In  terms  of  the  position  of  the  node  line  of  the  elastic 
mode.  When  z  ■>  0,  ©  -  0,  the  vertical  deflection  of  the  node  line  Is, 
by  definition,  zero.  Therefore  ${x,  y)  *  x  -  x^  over  the  lifting  sur¬ 
face  and  the  vertical  deflection  of  any  point  on  the  lifting  surface  Is 
z  +  x©  +  (x  -  x^)? 

while  the  slope  of  the  lifting  surface  Is  ©  +  £  . 


The  total  lift,  L,  and  the  moment  about  the  aerodynamic  center  of 
the  lifting  surface,  M,  car  be  expressed  In  terms  of  p(x,  y) j 
l  “  /p(x,  y)dS 

M  «  /  (x  -  xgc)p(x,  y)dS  -  /xp(x,  y)dS  -  xac  L 

It  follows  that  the  equations  of  motion  can  be  written 
mz  =■  L 
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where  x  »  x _  -  x  ,  .  The  lift  and  moment  on  the  lifting  surface  can 

an  ac  nt  % 

be  expressed  In  terms  of  the  vertical  displacement  of  the  aerodynamic 
center,  the  pitch  angle  of  the  surface,  and  the  derivatives  of  these 


quantities? 


za  c6. 
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In  matrix  forir,  the  eauatlons  of  motion  are 
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Introduce  c  as  a  unit  of  length,  ^  as  a  unit  of  Hire,  and  PSc  as  a 
unit  of  mass.  Also,  Introduce  the  radii  of  gyration 
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Then  the  following  are  convenient  non-dlrrenstona I  quantities! 
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The  equations  of  motion  can  now  be  written 


(IV-5) 

The  characteristic  equation  Is  obtained  from  the  determinant  of  the 
matrix  of  coefficients  In  (IV-5)! 


+  Aj  b5  +  A2  S2  +  A,  T  +  Aq)  -  0  ( IV-6) 
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where  u  and  «  are  defined  as  In  Section  III  of  this  report* 
oo 


In  order  to  determine  the  velocity  at  which  the  airframe  becomes 
neutrally  stable,  substitute  s  -  IQ  In  the  characteristic  equation, 
then  If  3  /  0» 

-  Aj  Q2  4  A,  -  0 
3^  -  Ag  a2  +  Ao  -  0 


or 


“  MsrJ+  Ao  "  0 


(IV-7) 


2  w 
This  equation  Is  linear  In  (-S'  If  we  divide  through  by  I  — 

[  O' 

u  la  given,  then  u  can  be  determined  from  the  ratio  - —  ,  Furthermore, 

G  O  <l> 

o 

since 


xac  CL 


P  S 


mr 


Is  a  linear  function  of  V  for  a  given  airframe.  It  follows  that  the 
velocity  at  neutral  stability  can  be  determined  from  (IV-7).  Note  that 
the  frequency  <i)Q  does  not  represent  the  short  period  frequency  of  the 
airframe  except  In  the  trivial  case  where  the  airframe  Is  rigid.  i»o 
Is  simply  a  convenient  way  of  specifying  the  velocity  V.  From  (IV-7) 
we  get 
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( JV-8) 


where 


(XV-9) 
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(IV-II) 


Equations  (IV-8),  (IV-9),  (IV-IO),  and  (XV-ll)  were  used  to  construct 

the  stability  boundaries  presented  In  the  Appendix.  The  frequency  of  the 

_  A. 

neutrally  stable  oscillation  Is  obtained  from  co  -  j 
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(IV-I2) 
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Six  parameters  suffice  to  determine  the  frequency  ratios  —  and 

o 

They  are 
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where  m  ts  the  mass  parameter,  m  ■  m/PSc.  These  parameters  are  the  same 
as  appeared  In  connection  with  the  simple  airframe  discussed  In  Section 


III  of  the  report,  except  that  the  parameter  — -  replaces  the  frequency 
ratio  parameter  — -  ,  and  the  node  line  position  xng  replaces  the  elastic 


“6 


axis  posl tlon  x 


In  the  limiting  case,  where  the  generalized  mess  of  the  elastic  mode 
Is  allowed  to  approach  zero  (F  »  0),  equation  (IV-6)  becomes  Identical 
with  equation  (III-I7)  of  Section  III,  provided  we  Identify  the  node  line 
position  with  the  elastic  axis  position  ><e  . 


As  was  found  previously,  the  presence  of  aerodynamic  damping  can 
have  a  profound  Influence  on  the  stability  of  the  airframe.  Consider, 
for  example,  the  case  where 


then 


Xan  ^ ~  * 


0t 


(Xv-13) 
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which  shows  that-  the  airframe  Is  neutrally  stable  at  zero  airspeed.  This 
Instability  Is  caused  by  the  loss  of  aerodynamic  damping  on  the  elastic 
mode  due  to  the  destabilizing  effect  of  C,  •  Thus,  since  structural  damp- 

q 

Ing  Is  Ignored,  the  structural  mode  Is  neutrally  stable  In  the  absence  of 
aerodynamic  damping.  Now  the  aerodynamic  damping  can  be  positive  or  nega¬ 
tive  In  general,  and  the  condition  (IV-13)  Is  just  the  condition  that  the 
aerodynamic  damping  be  zero.  If  -  0,  the  Instability  cannot  occur 

q 

since  equation  ( IV—  I  3)  cannot  be  satisfied.  This  type  of  Instability 
Is  entirely  similar  to  one  of  the  Instabilities  discussed  in  Section  III 

and  accounts  for  a  prominent  feature  of  the  stability  boundaries  present¬ 
ed  * 

ed  In  the  Appendix.  The  curves  of  .  given  In  the  Appendix  for  the 

“ofc 

case  of  subsonic  aerodynamics  display  a  sharp  dip.  The  value  of  xfln  at 
which  the  dip  occurs  Is  correctly  predicted  by  equation  (IV-13),  The 
curves  drawn  for  supersonic  aerodynamics  do  not  exhibit  this  feature, 
since  “0  for  these  curves.  Instability  due  to  loss  of  aerodynamtc 

damping  is  very  sensitive  to  the  presence  of  structural  damping  and  may 
be  masked  at  low  airspeeds  by  the  structural  damping, 

A  much  more  violent  Instability  will  be  exhibited  by  the  airframe 
of  Figure  lU  If  the  aerodynamic  center  Is  forward  of  the  node  line.  In 
this  case  the  short  period  frequency  and  the  elastic  mode  frequency  will 
come  close  together  at  a  sufficiently  large  airspeed.  At  higher  air¬ 
speeds,  one  of  these  two  roots  will  become  quite  unstable.  It  Is  useful 
to  examine  this  type  of  Instability  by  the  frequency  coalescence  technique. 

wofc  *s  *n  ®quati°n  (iv-19). 
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Ignore  all  damping  terms  In  the  equations  of  motion  (iV-i;).  Then  the 
characteristic  equation  can  be  obtained  from  (IV-6)  by  setting  C,  -  Cm 
■  0  and  dropping  the  terms  of  odd  order  In  5.  In  this  way  we  obtain  the 
following  characteristic  equation! 

S^+AS2*  B«0  (IV-II4.) 
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(iv- 1 5) 


(IV- 1 6) 


where  R  -  £-*•  , 
re 

o  2  2 

5  will  be  real  and  negative  unless  A  -  Is  negative.  If  A  —  i+B 

Is  negative  then  two  of  the  roots  of  equation  (lV-ll|)  will  have  a  positive 

real  part.  The  airframe  Is  therefore  on  the  verge  of  Instability  when 


A2  -  UB  •  0 


(IV- 17) 


or,  according  to  equations  (IV-15)  and  (IV-16),  when 

o  « 
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ac 
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(IV- IS) 


To  distinguish  the  value  of  ti)Q  determined  by  equation  (IV-8)  from 
the  value  given  by  equation  (IV-I8),  denote  the  solution  of  (TV-18)  by 


of  c 


Then 


Me  -1  /  „  xan 

-  “  I  +  V/  —  R - 

“ofc  '  xac 


(iv- 1 9) 


Thus,  according  to  the  frequency  coalescence  criterion,  there  are  only 
two  parameters  Involved  In  the  stability  problem. 
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The  parameter 


R 


( general  I  zed  mass  ratio) 


xan  B  xac  "  xni 
xac  xao 

Is  negative  If  the  aerodynamic  center  Is  forward  of 

the  node  1 1 ne.  According  to  frequency  coalescence.  Instability  can  occur 
only  If  the  aerodynamic  center  Is  forward  of  the  node  line  of  the  elastic 
wofc  xan 

mode.  -  Is  plotted  versus  -  In  Figure  15  for*  several  values  of 

ue  xac 

R.  The  frequency  coalescence  solution  of  the  stability  problem  Is  re¬ 
markably  simple,  however  It  can  be  highly  unconservative  due  to  the  fact 
that  the  effects  of  aerodynamic  damping  are  entirely  Ignored.  Also,  only 
one  of  the  two  essenttally  distinct  Instability  mechanisms  Is  explained 
by  the  frequency  coalescence  approach. 


The  discrepancies  between  the  predictions  of  the  frequency  coales¬ 
cence  method  and  predictions  of  the  more  exact  theory  are  clearly  exhibit¬ 
ed  In  the  stability  boundaries  of  the  Appendix.  On  the  graphs  presented 
In  the  Appendix  the  frequency  coalescence  solution  Is  simply  the  straight 
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Figure  I Flexible  Airframe  With  a  Rigid  Lifting  Surface 


Figure  15.  Stability  Boundaries  Cons 


V.  CLASSIFICATION  OF  AEROELASTIC  STABILITY  PROBLEMS 


The  study  and  analysts  of  the  stability  of  airborne  vehicles  was 
historically  divided  Into  three  fields  of  Interest: 

1.  Stability  and  control. 

2.  Flutter. 

3.  Steady  state  aeroelastlctty  (divergence). 

A  definition  of  scopes  of  these  fields  has  never  been  established;  how¬ 
ever,  In  their  origin  these  areas  of  technology  were  directed  at  specific 
phenomena. 

1.  The  analytical  effort,  most  commonly  called  the  study  of 
"stability  and  control",  was  concentrated  on  aircraft 
stability  problems  which  exist  when  the  aircraft  Is  a  rigid 
body.  The  effects  of  elasticity  have  been  Incorporated  Into 
these  analyses  but  often  In  such  a  manner  as  to  exclude  the 
dynamic  response  of  the  system  In  any  of  Its  finite  frequency 
normal  modes. 

2.  The  field  of  flutter  analysis  has  been  limited  to  the  study  of 
aeroelastlc  Instabilities  which  arise  from  the  "coupling"  of 
two  or  more  finite  frequency  elastic  modes  of  the  aircraft. 

In  many,  and  possibly  most,  flutter  analyses  of  aircraft,  the 
vehicle  has  been  represented  by  a  relatively  small  number  of 
(10  or  less)  finite  frequency  elastic  modes  of  the  free  system. 
Zero  frequency  modes  have  commonly  been  omitted  from  flutter 
analyses  along  with  consideration  of  the  generalized  flexibility 
of  modes  higher  than  those  specifically  Included  In  the  analyses. 

3.  Steady  state  aeroelastlc  Instability,  commonly  called  divergence. 
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has  been  a  recognized  problem  since  the  very  early  days  of 
aircraft  design.  The  analyses  of  steady  state  aeroelasttc 
phenomena  In  essence  Ignore  the  Inertia  of  the  vehicle  except 
for  the  consideration  of  "Inertia  relief"  In  free  systems. 

Recent  years  have  brought  about  some  changes  In  the  analytical  fields 
discussed  above  which  tend  to  merge  their  technologies.  The  results  of 
this  study  Indicate  that  a  unification  of  these  analytical  fields  Is 
desirable  and  necessary  to  produce  a  reliable  aeroelastlc  stability 
analysts  of  a  modern  airborne  system. 

CLASSIFICATION  OF  AEROELASTIC  PHENOMENA 

Collar  presented  a  classification  of  aeroelastlc  problems  In  the 
well  known  "Collar’s  triangle  of  forces"  by  which  aeroelastlc  phenomena 
are  described  as  being  the  result  of  various  combinations  of  aerodynamic, 
elastic  and  Inertia  forces.  We  propose  another  aeroelastlc  triangle 
(Figure  16)  for  free  systems  as  a  supplement  to  Collar's  triangle. 

A  classification  of  aeroelastlc  phenomena  Is  made  In  the  triangle 
of  Figure  |6  by  the  normal  modes  of  the  system  which  are  required  In  a 
representation  of  the  system  to  produce  the  phenomenon.  The  apexes  of 
the  triangle  of  Figure  16  aret 

0  -  zero  frequency  or  rigid  body  modes 
F  -  finite  frequency  modes 

oo  -  Infinite  frequency  modes  or  residual  flexibility 
Since  Figure  16  applies  to  aeroelastlc  problems,  aerodynamic  forces  are 
Involved  In  all  phenomena  described. 

The  problems  Included  In  the  nunbered  boxes  are* 
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1.  Classical  stability  and  control  problems  Involving  a  rigid 
vehicle  and  aerodynamic  forces. 

2.  Classical  flutter  problems  Involving  2  or  more  finite  frequency 
modes  and  aerodynamic  forces. 

3.  Classical  steady  state  aeroelastlc  problems  Involving  only 
elasticity  and  aerodynamic  forces. 

U.  "Mode  Interaction.”  To  a  stability  and  control  analyst  this 
problem  Is  described  as  one  where  elastic  modes  have  a  pro¬ 
nounced  effect  on  the  low  frequency  response  of  the  system. 

To  a  flutter  analyst  this  Is  a  problem  Involving  coupling  of 
the  rigid  body  modes  with  the  elastic  modes.  This  phenomenon 
Is  shown  to  be  a  stability  problem  In  the  analysis  of  Section 
IV  of  this  report  and  Is  undoubtedly  a  problem  In  the  predic¬ 
tion  of  the  dynamic  behavior  of  a  system  at  a  velocity  below 
the  flutter  speed. 

The  analysis  of  Section  IV  shows  that  this  stability 
problem  can  be  predicted,  for  many  systems,  by  "frequency 
coalescence”  methods  (omitting  damping  terms)  but  for  certain 
systems  the  damping  terms  can  have  a  dominating  effect  on  the 
so  I utlon. 

5.  "Elastic  Interaction."  This  problem  area  will  Include  the 

modification  of  "rigid  body"  dynamic  response  by  the  elasticity 
of  the  system  and  the  stability  problem  shown  by  the  analysis 
of  Section  III  of  this  report  where  a  free  system  can  be  un¬ 
stable  by  virtue  of  only  plunge  or  pitch  flexibility  of  its 
aerodynamic  surface. 
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This  stability  problem  Involving  plunge  flexibility  Is 
quite  Interesting  because i 

a.  It  Is  a  flutter  problem  which  can  exist  for  on 
Idealized  system  which  has  no  finite  frequency 
normal  modes. 

b.  It  Is  a  flutter  problem  for  which  the  ’’mechanism” 
of  Instability  can  be  easily  understood. 

Most  flutter  phenomena  are  the  result  of  complicated,  relative 
phasing  between  structural  motions  which  cause  aerodynamic 
forces  to  either  add  or  subtract  energy  from  the  system.  For 
this  case  of  "elastic  Interaction"  the  structural  deflections 
have  no  relative  phase  angle  and  the  phenomenon  can  be  easily 
explained  as  follows. 

Consider  the  airframe  of  Figure  13  with  a  rigid  pitch 
spring  (Kg)  and  a  flexible  plunge  spring  (K^).  Figure  17 
shows  the  airframe  at  the  Instant  of  maximum  deflection  (of 
unit  pitch  angle)  Jn  Its  short  period  mode  for  various  values 
of  velocity. 

Figure  17a  corresponds  to  a  low  value  of  velocity.  The 
short  period  mode  frequency  Is  low  as  well  as  the  lift  force 
L.  The  short  period  mode  Is  well  damped  by  virtue  of  the 
lift  due  to  the  plunge  velocity  h. 

Figure  17b  corresponds  to  a  higher  value  of  velocity. 

The  lift  force  has  Increased  from  that  of  Figure  I7aj  and, 
therefore,  the  short  period  mode  frequency  and  the  deflection 
of  the  spring  have  Increased.  Since  h  Is  smaller,  the 
damping  of  the  short  period  mode  has  decreased. 
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Figure  l?c  corresponds  to  a  higher  velocity  than  Figure 
17b*  At  this  speed  the  lift  force  on  the  aerodynamic  surface 
Is  sufficient  to  deflect  the  plunge  spring  («h)  the  distance 
y  and  therefore  reduce  the  plunge  of  the  surface  to  zero. 
Since  the  plunge  of  the  surface  Is  zero,  the  damping  due  to 
plunge  Is  also  zero.  If  this  system  Is  subjected  to  a  still 
higher  velocity  the  aerodynamic  surface  will  plunge  In  a 
direction  opposite  from  Its  attachment  point.  Since  the 
plunge  velocity  h  of  Figure  17  contributed  positive  damping 
to  the  system,  It  Is  only  reasonable  that  when  the  algebraic 
sign  of  h  Is  changed  the  damping  will  be  negative.  Since  the 
damping  In  the  short  period  mode  Is  attributable  primarily  to 
the  plunge  velocity,  the  system  will  be  unstable  at  a  dynamic 
pressure  slightly  higher  than  that  depicted  In  Figure  l/c. 

6.  "General  Aeroelastlc  Problem."  Section  II  of  this  report 
presented  an  analysis  and  comparison  of  the  solution  of  con¬ 
figuration  1*  of  reference  I.  It  has  been  demonstrated  that 
a  reliable  prediction  of  the  dynamic  response  of  this  con¬ 
figuration  can  be  attained  only  through  the  Inclusion  of  zeno 
frequency  modes,  finite  frequency  modes,  and  the  residual 
flexibility  of  all  higher  modes  of  the  system  In  the  modal 
representation. 

The  analysis  of  Section  II  omitted  the  aerodynamic  damp¬ 
ing  terms  for  simplicity.  This  assumption  was  apparently 
reasonable  for  the  solution  of  configuration  4  because  of  the 
relatively  high  altitude  considered.  The  analysis  of  Section 
III  shows  that  this  assumption  cannot  be  used  In  the  general 
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case  and  for  many  configurations  can  lead  to  an  extremely 
unconservative  prediction  of  the  flutter  speed. 

7«  The  problem  area  described  by  area  7  of  Figure  16  can  be 
described  as  an  aeroelastlc  Instability  which  can  be  repre¬ 
sented  analytically  when  only  an  elastic  mode  and  the  residual 
flexibility  of  the  system  are  Included.  Such  a  phenomenon  ts 
not  known  to  exist  but  It  Is  postulated  that  the  problem  area 
will  be  discovered  In  the  future. 
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Figure  1 6 


I _ I 


.  Diagram  of  Aeroelastlc  Phenomenon 
of  Free  Systems 
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VI.  COMPUTER  STUDIES  OF  MISSILE  CONFIGURATION 
AND  COMPARISON  WITH  THEORY 


Studies  were  made  on  the  Computer  Engineering  Associates'  passive 
analog  computer  of  a  number  of  "airborne  missile"  configurations.  The 
purpose  In  making  these  computer  studies  was  to  provide  a  check  or 
Independent  verification  of  the  analyses  of  this  report  and  to  study 
the  effects  of  parameters  not  Included  In  these  analyses. 

The  basic  missile  configuration  considered  In  all  computer  studies 
was  configuration  2  of  reference  I.  The  geometry,  mass  and  stiffness 
data  for  this  basic  configuration  are  repeated  In  Figures  18,  19,  20 
and  Table  I.  During  the  conduct  of  this  study  numerous  variations  were 
made  In  this  conf Iguratlon.  In  most  computer  analyses  the  forward  aero¬ 
dynamic  surface  was  omitted.  All  variations  to  the  configuration  are 
summarized  In  Table  2.  In  all  cases  the  aerodynamic  surfaces  were 
considered  to  be  rigid  and  elastically  restrained  by  a  "pitch"  and 
"plunge"  spring  at  the  missile  station  listed  In  Table  2  as  "e.a. 
Station". 

The  aerodynamic  forces  were  represented  as  fol lowst 

1,  The  forces  on  the  body  of  the  missile  were  assumed  to 
be  zero. 

2.  The  lift  forces  on  the  surfaces,  normal  to  the  surfaces 
and  positive  upward  are  given  by 


where  Z  Is  the  normal  deflection  (positive  up)  and  6  Is 
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the  pitching  slope  (positive  nose  up)  of  the  surface  at 
the  aerodynamic  center  of  the  surface. 

3.  The  moments  about  the  centers  of  pressure  of  the  surface 
are  alven  by« 

g 

U  ^  P  V  S  ^  s  0 
i+.  The  effects  of  downwash  were  neglected. 

The  "baste"  flight  condition  Is  defined  by  the  following  aerodynamic 
constants* 

Velocity  *  2250  mph 

Al ti tude  ■  sea  level 

Dynamic  Pressure  ■  qQ  »  90  lb./ln.2 

Lift  Curve  Slope  •  *»  1.5  per  rad« 

a 

The  flutter  speeds  measured  In  the  studies  were  tabulated  as  a  fraction 
of  the  basic  velocity  (Vq),  thus  flutter  speed  Is  listed  as  the  ratio 
vAo. 

The  computer  studies  of  the  configurations  listed  In  Table  2  deter¬ 
mined  the  dynamic  pressure  corresponding  to  flutter.  In  all  cases  the 
flutter  encountered  was  of  the  general  type  discussed  In  this  report* 
that  Is,  It  Involved  the  rigid  body  modes  of  the  system  as  Indicated  by 
a  flutter  frequency  below  the  lowest  natural  frequency  of  structure  In 
a  vacuum. 

RIGID  FUSELAGE  -  CASES  I  THROUGH  9 

Cases  I  through  9  considered  a  rigid  fuselage  with  a  single  aero- 
dynamte  surface  elastically  attached  near  the  aft  end  of  the  missile. 
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Listed  at  the  right  side  of  Table  2  are  the  values  of  flutter  speed 
measured  on  the  analog  computer,  calculated  from  the  crtterla  of  Section 
III  and  from  the  criteria  of  Section  IV, 

Cases  I  through  5  considered  an  aerodynamic  center  which  coincided 
with  the  elastic  axis  of  the  lifting  surface. 

Cases  I,  2  and  5  are  identical  except  for  the  pitch  Inertia  of  the 
elastically  supported  aft  surface.  The  correlation  between  theory  and 
measurement  was  fair  for  cases  I  and  5  where  the  Inertia  was  quite  small 
(zero  In  case  5)  hut  much  poorer  In  case  2  where  the  Inertia  was  10 
times  larger.  The  flutter  speed  calculated  from  the  criteria  of  Section 
III  Is  Identical  for  cases  I,  2  and  5  because  the  mass  of  the  surface 
Is  completely  Ignored  In  the  analysis  of  Section  III.  The  discrepancy 
between  the  theory  of  Section  III  and  measurement  In  case  2  Is  easily 
explained  by  omission  of  the  mass  from  the  Section  III  criteria  but 
the  discrepancy  between  the  measurement  and  the  criteria  of  Section  IV 
cannot  be  explained. 

Cases  3  and  1+  are  Identical  except  that  the  mass  of  the  lifting 
surface,  elastically  mounted  by  a  plunge  spring.,  was  omitted  In  case  it. 
The  Identical  flutter  speed  was  measured  on  the  computer  for  these  2 
cases  as  also  was,  of  course,  predicted  by  the  theory  of  Section  III 
which  Ignored  the  mass  completely.  The  crtterla  of  Section  IV  could 
not  be  applied  to  cases  it  and  5  due  to  the  absence  of  a  finite  frequency 
elastic  mode.  The  correlation  between  theory  and  measurement  Is  con¬ 
sidered  reasonably  good  for  these  cases.  The  discrepancies  can  be 
ascribed  to  the  equivalent  structural  damping  inherently  present  In 
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the  passive  analog  computer. 


Cases  k  find  5  are  Interesting  specimens  of  flutter  because  In  a 
vacuum  the  structures  In  both  cases  have  no  finite  frequency  modes. 

When  described  by  their  normal  coordinates  these  configurations  possess 
two  zero  frequency  degrees  of  freedom  and  one  Infinite  frequency  degree 
of  freedom.  Nevertheless  these  configurations  exhibit  a  flutter  In¬ 
stability  at  a  flntte  frequency. 

Cases  6  and  7  ere  similar  to  cases  I  and  3 »  respectively,  except 
that  the  aft  aerodynamic  surface  was  shifted  forward  10#  of  Its  chord 
In  cases  6  and  7  while  retaining  the  attachment  point  of  the  surface 
to  the  fuselage  at  the  same  point  on  the  fuselage.  The  effective 
elastic  axis  of  the  aerodynamic  surface  was,  therefore,  at  the  60# 
chord,  10#  of  the  chord  aft  of  the  aerodynamic  center.  The  configura¬ 
tions  of  case  6,  where  the  pitch  restraint  between  the  fuselage  and  the 
aerodynamic  surface  Is  flexible,  tends  toward  steady  state  elastic 
divergence.  Divergence  would  be  encountered  at  a  dynamic  pressure 
corresponding  to  V/VQ  of  1.39.  Table  2  shows  that  a  dynamic  Instability 
Is  encountered  at  a  somewhat  lower  speed  than  the  divergence  speed. 

The  comparison  of  cases  6  and  7  to  cases  I  and  3  demonstrate  the 
sensitivity  of  the  flutter  phenomenon  to  the  relative  locations  of  the 
aerodynamic  center  and  the  elastic  axis  of  the  lifting  surface.  This 
extreme  sensitivity  may  be  considered  as  an  explanation  for  some  of  the 
discrepancies  In  theoretical  predictions  and  computer  measurements  shown 
In  Table  2.  The  precise  location  of  the  aerodynamic  center  Is  difficult 
to  control  In  the  experimental  procedure  followed  on  the  analog  computer. 
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It  seems  also  worthwhile  to  point  out  that  the  precise  location  of  the 
aerodynamic  center  Is  generally  not  known  for  a  physical  system  and.  In 
fact,  cannot  be  defined  as  physical  point  except  for  Idealized  systems. 

Cases  8  and  9  are  similar  to  cases  7  and  6,  respectively,  except 
that  the  lifting  surface  was  moved  aft  J00  Inches  In  cases  8  and  9. 

This  variation  affected  the  flutter  speed  very  little  In  the  computer 
measurements  and  had  a  somewhat  larger  effect  on  the  theoretical  pre¬ 
dictions. 

FUSELAGE  WITH  SINGLE  BENDING  DEGREE  OF  FREEDOM  -  CASE  10 

Case  10  considered  a  configuration  where  the  entire  system  was 
rigid  except  for  a  single  flexibly  restrained  "hinge"  at  fuselage  sta¬ 
tion  700.  The  purpose  of  studying  thts  design  was  to  experimentally 
evaluate  the  phenomenon  shown  on  the  charts  of  Section  IV  where  the 
flutter  speed  Is  zero  for  specific  parameter  combinations.  The  flexi¬ 
bility  of  the  system  was  reduced  to  the  single  flexibly  restrained 
"hinge"  In  the  fuselage  to  preclude  considerations  of  "residual  flexi¬ 
bility".  The  location  of  the  lifting  surface  was  varied  In  thts  case 

x  .  -  x 

to  provide  a  variation  In  the  dimensionless  parameter  — — - —  , 

ac 

Figure  21  shows  a  comparison  between  the  flutter  speeds  predicted 
by  the  theory  of  Section  IV,  frequency  coalescence,  and  analog  computer 
measurements.  The  radical  difference  between  theoretical  prediction  by 
Section  IV  and  computer  measurements  Is  easily  explained  by  the  presence 
of  a  small  amount  of  equivalent  structural  damping  In  the  computer 
sol utlon. 
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Figure  22  shows  some  typical  V-g  diagrams  for  the  configuration  of 
case  10.  The  curve  labeled  "structural  damping"  shows  the  equivalent 
structural  damping  level  which  existed  In  the  analog  computer  at  the 
values  of  dimensionless  velocity  v//^.  For  the  curve  corresponding  to 


x  .  -  x 
ni  ac 


ac 


of  3.8lf  the  Intersection  with  the  structural  damping  curve 


Is  at  a  value  of  V/A/q  virtually  Identical  to  that  at  Its  Intersection 
with  the  0  damping  axis  which  Is  considered  to  be  the  flutter  speed. 
*nl  "  xac 

At  this  value  of  —  ■—  theory  and  computer  measurement  are  In 

xnl 

reasonable  agreement  as  shown  In  Figure  21. 


nl  oe 

The  curve  on  Figure  22,  corresponding  to  ■■■  . . .  ■■  of  .261+  never 

xae 

Intersects  with  the  "structural  damping"  curve*  and  If  the  structural 
damping  curve  were  subtracted  out  of  the  solutions,  this  configuration 
would  apparently  be  unstable  at  ell  finite  velocities  as  shown  In 
Figure  21. 


This  single  example  cannot  be  considered  to  prove  that  the  low 
speed  Instability  problem  shown  in  Section  IV  can  be  completely  Ignored 
In  the  presence  of  a  small  amount  of  structural  damping.  This  example 
does  show  that  this  problem  may  not  be  as  serious  as  might  be  concluded 
from  a  consideration  of  only  the  results  of  Section  IV. 


EFFECT  OF  ALTITUDE  -  CASES  II  AND  12 

Cases  II  and  12  may  be  compared  with  cases  I  and  3,  respectively, 
to  show  the  effect  of  an  Increase  In  altitude.  Cases  II  and  12  correspond 
to  an  air  density  l/lO  of  the  air  density  considered  In  cases  I  and  3* 


70 


The  change  In  the  flutter  speed  Is  shown  to  be  an  Increase  of  12$  to 
16$,  It  may  be  concluded  that  the  "mode  Interaction"  phenomenon  Is 
not  sensitive  to  air  density  or  altitude  parameters. 

FLEXIBLE  FUSELASE  -  CASES  13  THROUGH  25 

The  comparison  between  theory  and  computer  measurements  for  cases 
Involving  a  flexible  fuselage.  In  all  cases,  showed  rather  poor  cor¬ 
relation.  This  lack  of  correlation  Is  probably  due  tot 

a.  The  presence  of  higher  elastic  modes  In  the  computer  repre¬ 
sentation  which  were  Ignored  In  the  theoretical  analyses. 

b.  Equivalent  structural  damping  In  the  computer  measurements 
which  was  Ignored  In  the  theoretical  analyses. 

The  analyses  of  Sections  II  and  III  were  not  applied  to  these  con¬ 
figurations  because  they  were  known  to  Include  assumptions  which  pre¬ 
cluded  their  applicability.  The  analysis  of  Section  II  cannot  be  con¬ 
sidered  applicable  to  any  cases  of  this  configuration  because  It  Ignores 
aerodynamic  damping.  The  results  of  cases  I  through  9  and  the  analysis 
of  Section  III  Identify  the  mode  Interaction  phenomenon  of  this  "missile 
configuration"  to  be  Intimately  related  to  the  aerodynamic  damping.  How¬ 
ever  the  flutter  phenomenon  observed  In  some  cases  of  this  configuration 
were  undoubtedly  related  to  the  Instability  studied  In  Section  II. 

EFFECT  OF  FLEX  I BLV  MOUNTED  LIFTING  SURFACE  -  CASES  13  THROUGH  16 

Cases  13  and  12;  consider  configurations  that  are  Identical  except 
for  the  flexibility  of  the  attachment  of  the  lifting  surface  to  the 
fuselage.  Case  1 3  considered  the  attachment  to  be  rigid  and  case  ll;  con¬ 
sidered  the  attachment  In  pitch  and  plunge  to  be  of  "basic"  flexibility. 
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The  flutter  speeds  measured  for  these  cases  were  virtually  Identical, 
but  the  system  damping  at  speeds  less  than  the  flutter  speed  were  quite 
different.  Figures  23  and  2k  show  the  V-g  diagrams  for  these  cases  In¬ 
cluding  the  "short  period  mode"  and  the  lowest  "elastic  mode  of  the 
system".  In  case  13,  where  the  aerodynamic  surface  was  rlgtdly  connect¬ 
ed  to  the  fuselage,  the  short  period  mode  remained  well  damped  throughout 
the  speed  range  up  to  the  vicinity  of  the  flutter  speed  where  Its  damping 
Increased.  In  case  ll;,  where  the  lifting  surface  was  elastically  connect¬ 
ed  to  the  fuselage,  the  damping  In  the  short  period  mode  became  progres¬ 
sively  smaller  as  velocity  Increased  until  It  became  unstable  at  the 
f I  utter  speed. 

The  simplified  analyses  of  this  study  are  not  capable  of  either 
predicting  or  providing  an  explanation  for  this  difference  between  cases 
13  and  lU.  The  results  of  this  study  are  capable,  however,  of  defining 
the  minimum  complexity  of  an  analysis  which  would  be  applicable  to  these 
cases  as  an  analysis  which  Includes  consideration  of 

a.  Two  rigid  body  modes, 

b.  One  elastic  mode, 

c.  The  residual  f lex! bl 1 1 ty  of  the  higher  modes  of  the  system. 

d.  Aerodynamic  damping  terms. 

e.  Prediction  of  system  damping  at  speeds  less  than  the  flutter 
speed. 

Such  an  analysis  of  a  specific  configuration  is  not  considered  overly 
complex  for  practical  solution,  but  the  complexity  does  appear  too  great 
for  the  type  of  generally  applicable  solutions  sought  tn  this  study. 

It  should  be  noted  that  the  Intention  of  the  analyses  of  this  study 


72 


was  to  predict  the  susceptibility  of  a  configuration  to  "mode  Inter¬ 
action"  and  that  the  difficulties  experienced  In  predicting  precise 
flutter  speeds  or  the  response  of  the  system  at  subcritlcal  speeds 
does  not  Impair  the  fulfillment  of  this  Intention. 

A  comparison  of  cases  1 5  and  1 6  show  that  when  the  rigidly  mounted 
surface  Is  moved  aft  of  the  node  line  of  the  first  elastic  mode,  the 
Instability  disappears}  but  that  when  the  flexible  mounting  of  the  lift¬ 
ing  surface  Is  Introduced  the  system  ts  again  unstable  at  virtually  the 
same  speed  as  In  cases  1 3  and  1 4  where  the  surface  Is  forward  of  the 
node  line. 

EFFECTS  OF  CANARD  SURFACE  -  CASES  18  AND  19 

The  configurations  studied  In  cases  18  and  I?  are  Identical  to 
those  considered  In  cases  14  and  13.  respectively,  except  that  the 
canard  surface  of  the  "basic"  configuration  Is  added  In  cases  10  and 
19.  The  effects  of  the  canard  surface  are  simply  to  Increase  the 
flutter  speed  about  7%,  The  V-g  curves  of  cases  18  and  19  are  similar 
In  character  to  those  of  cases  l4  and  13.  respectively,  shown  In  Figures 
23  and  24, 

EFFECTS  OF  THE  POSITION  OF  THE  SINGLE  LIFTING  SURFACE  -  CASES  22  THROUGH  23 

Cases  13  through  |6  and  cases  22  through  25  show  the  effect  of  posi¬ 
tioning  the  lifting  surface  at  four  locations  along  the  fuselage  from 
station  800  to  station  1000,  At  each  location  the  effect  of  flexibility 
of  the  lifting  surface  attachment  was  Investigated.  As  shown  In  Table 
2,  the  flutter  speed  was  affected  very  little  by  the  variations  in  position 
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of  the  lifting  surface  or  the  flexibility  of  the  aerodynamic  surface 
attachment.  In  each  case,  when  the  attachment  was  rigid  the  short 
period  mode  damping  remained  constant  over  the  speed  range  studted, 
and  when  the  attachment  was  flexible  the  short  period  mode  damping 
became  progressively  small  as  speed  was  Increased  until  the  short 
period  mode  became  unstable  at  the  flutter  speed. 
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figure  19.  Body  Flexibility  of  Configuration  2 
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Figure  21.  Comparison  of  Computer  Results  and  Theory  -  Case  10 


figure  22.  Effect  of  Parasitic  Computer  Damping  Upon 
Neutral  Stability  Determination  -  Case  10 
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Figure  23.  V-g  Curve  for  Case  1 3 
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FI  gup*  2lu  V-g  Curve  for  Case  lit 
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Table  I.  Mass  Distribution  for  Configuration  2 


Mass  No. 
(Figure  20) 

Station 

Condition  3° 
Mass 

Condition  31 

Mass 

(In.) 

1 

0 

4.728 

4.728 

2 

50 

18.911 

18.911 

3 

150 

18.911 

18.911 

4 

250 

31.605 

31.605 

5 

350 

56.734 

18.911 

6 

U50 

82.122 

27.374 

7 

550 

94.616 

31.605 

8 

650 

94.816 

31.605 

9 

750 

51.812 

17.270 

10 

900 

51. 812 

17.270 

it 

1000 

189.632 

63.210 

12 

MOO 

45.854 

15.285 

13 

1200 

22.929 

7.642 

Forward 

Surface 

250 

1.554 

1.554 

Aft 

Surface 

900 

15.544 

15.544 

PI tchlng  Mass  Moment  of  Inertia 

■ 

2 

lb. -sec.  -In. 

2 

fb.-sec.  -In* 

Forward 

Surface 

250 

155.4 

155.4 

Aft 

Surface 

900 

2487. 

2487. 
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Table  2, 


Table  2.  Missile  Configuration  -  Description  of  Cases 


Aft  Surface 

Air 

Density 

q/{^  at  Flutter 

e.a.  Position 

In  #  Chord 

a.c.  Position 
In  $  Chord 

Area 

Plunge 

1 nertla 

Pitch 

1  nertla 

Plunge 

P  lex. 

Pitch 

Flex. 

Measured 

CAlc.  Ill 

Calc.  IV 

50# 

50$ 

Basic 

Basic 

Sasic 

Rigid 

Basic 

Basic 

2.56 

2.80 

2.73 

■ 

■ 

Basic 

x  10 

Rigid 

Baste 

■ 

■ 

1.66 

2.80 

2.2l* 

Basic 

Basic 

Baste 

Rigid 

2.26 

2.13 

2.18 

0 

Basic 

Basic 

Rigid 

2.26 

2.13 

50$ 

Basic 

0 

Rigid 

Basic 

2.69 

2.80 

60$ 

■ 

■ 

Basic 

Rigid 

Basle 

l.ll 

1.19 

1.22 

i 

■ 

■ 

Baste 

Rigid 

2.20 

2.20 

2.21* 

1 

Basic 

Rigid 

2.20 

2.11 

1.86 

60$ 

50$ 

■ 

■ 

Rigid 

Basic 

1 . 10 

M3 

1.17 

25 $ 

25# 

Rigid 

Rigid 

Basic 

50$ 

50# 

_ 

Basic 

x  ,  1 

2.97 

■ 

Basic 

Rigid 

x  .  1 

2.53 

■ 

■ 

Rigid 

Basic 

1. 10 

.566 

■ 

m 

■ 

1.12 

■ 

■ 

i 

■ 

■ 

Basic 

Basic 

■ 

■ 

1.16 

50$ 

50# 

Basic 

x  .  l6l 

Rigid 

Rigid 

Basic 

Stable 

50$ 

50# 

Basic 

Baste 

x  .  l6l 

Bas  I  c 

Basic 

1. 18 

50$ 

50# 

Basic 

Basle 

x  .  I6l 

Basic 

1.17 

50$ 

50# 

Basic 

Basic 

x  .  l6l 

Rigid 

Rigid 

Baste 

1.08 

.546 

m 

■ 

■ 

1 

■ 

Basic 

Basic 

■ 

■ 

1.09 

.493 

m 

■ 

■ 

■ 

■ 

■ 

Basic 

Basic 

■ 

■ 

1.18 

.549 

50$ 

50# 

Baste 

Basic 

m 

Rigid 

Rigid 

Basic 

.920 
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VII.  CONCLUSIONS 


1.  An  accurate  approximation  for  the  aeroelastlc  behavior  of  a  structure 
may  be  made  In  terms  of  Its  normal  coordinates  by  Including  some  of 
Its  modes  explicitly  and  the  "residual  flexibility"  approximation  to 
all  higher  modes.  This  approximation  for  the  system,  derived  In 
reference  I,  Is  applicable  In  the  presence  of  "mode  Interaction", 

This  conclusion  Is  drawn  from  the  fact  that  "mode  Interaction"  was 
found  to  depend  on  only  the  parameters  listed  In  the  foregoing 
approximation. 

2.  In  the  presence  of  "mode  Interaction,"  the  accurate  determination 
of  "residual  flexibility"  Is  extremely  Important.  In  the  presence 
of  "mode  Interaction"  the  dynamics  of  a  system  can  be  very  sensitive 
to  slight  changes  In  flexibility. 

3.  Systems  which  tend  toward  steady  state  divergence  are  particularly 
susceptible  to  "mode  Interaction." 

i4«  In  most  free  systems  where  steady  state  divergence  Is  predicted  to 
occur,  a  mode  Interaction  or  flutter  Instability  will  probably  occur 
at  a  velocity  lower  than  that  predicted  for  steady  state  divergence. 

5.  Aerodynamic  damping  terms  can  have  a  destabilizing  effect  on  a  free 
system. 

6.  Aeroelastlc  systems  which  possess  no  finite  frequency  normal  modes, 
can  be  susceptible  to  a  finite  frequency  flutter  Instability, 


7.  A  reliably  accurate  method  tor  predicting  the  dynamic  response  of 
a  general,  free,  aeroelastlc  system  must  Include  consideration  of 


a.  The  rigid  body  modes  of  the  system, 

b.  The  lowest  elastic  mode  of  the  system. 

c.  The  "residual  flexibility"  of  all  higher  modes. 

d.  Aerodynamic  damping. 

e.  Structural  damping. 

8,  Flutter  analyses  of  free  systems  must  Include  the  zero  frequency 
modes  of  the  systems. 
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APPENDIX 


STABILITY  BOUNDARIES  FOR  AN  AIRFRAME  HAVING  TWO 
RIGID  BODY  MODES  AND  ONE  ELASTIC  MODE 

A  brief  summary  of  the  analysis  of  section  IV  Is  first  presented. 

The  aeroelastlc  system  of  Figure  25  Is  characterized  by  two  rigid 
body  modes  and  a  single  elastic  mode.  It  Is  assumed  that  aerodynamic 
forces  act  only  upon  a  rigid  lifting  surface  and  that  these  forces  produce 
a  lift  and  moment  about  the  aerodynamic  center  which  are  given  by 


where 

i  p 

q  »  ^  P  V  Is  the  dynamic  pressure. 

Is  the  lift  coefficient, 
a 

CL 

1  -  2x  •  jp- 1  Is  the  ratio  of  lift  coefficient  caused  by  pitch 
CL 

a 

rate  to  the  lift  coefficient  caused  by  angle  of  attack, 

C  Is  the  pitch  damping  moment  coefficient. 

"q 

c  Is  the  chord  length. 

s  Is  the  Laplace  transform  differentiation  operator  • 

The  equations  of  motion  for  the  system  of  Figure  25  can  be  written 

mz  o  L 
2  •• 

mr  0-M+x  L 

fiC 

m  re2(?+  Q  »  M  ♦  (xac  -  x^L 

and  the  motion  of  the  lifting  surface  Is  related  to  the  normal  coordinates 
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2,  6,  4  as  fol  lows: 


zac  “  z  +  xac  6  +  (xac  “  xn^ 


6ac  *  ®  *  « 


The  symbols  appearing  In  these  equations  have  the  following 
significance: 

m  is  the  total  mass  of  the  system  (generalized  mass  of  rigid  plunge 

mode) . 

2 

m  r  Is  the  total  pitching  Inertia  about  the  center  of  mass  (general¬ 
ized  mass  of  rigid  pitch  mode). 

p 

m  rfl  Is  the  generalized  mass  of  the  elastic  mode  when  the  elastic 
mode  shape  Is  normalized  by  making  the  slope  of  the  elastic  mode 
unity  at  the  node  line. 

p 

Is  the  frequency  of  the  elastic  mode, 
x  Is  the  x-coordlnate  of  the  aerodynamic  center  relative  to  the 
coordinate  system  of  Figure  25  (If  the  aerodynamic  center  Is  behind 
the  center  of  mass,  then  x  wl  1 1  be  negative). 

aC 

xn<  Is  the  x-coordlnate  of  the  node  line  of  the  elastic  mode  (negative 
If  the  node  line  Is  behind  the  center  of  mass). 


Analysis  of  the  above  system  of  equations  produces  the  following 
resul  ts: 


a.  Let  biQ  be  a  fictitious  short  period  frequency  defined  by 


(jPvS)SCL  *, 


ac 


PSCL  xac 
a 


m  r 


2  m  r 


Then  the  system  Is  neutrally  stable  whenever  the  frequency  ratio  — 

e 
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takes  on  the  value  given  by 


b.  As  a  corollary,  It  follows  that  the  stability  of  the  system  of 
Figure  25  Is  determined  by  six  parameters: 


nl 


ac 


ac 


R 


r 


x 


ac 


ac 
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The  last  of  these  parameters  Is  much  less  Important  than  the  other  five, 
and  the  first  two  parameters  are  decisive  when  the  effects  of  damping 
are  negligible. 

c.  If  all  damping  terms  are  Ignored  In  the  equations  of  motion, 
then  a  frequency  coalescence  type  of  analysis  can  be  applied,  with  the 
following  result.  Let  uq^c  be  the  value  of  to0  at  which  the  system  Is 
neutrally  stable  according  to  frequency  coalescence.  Then 

“ofc 
to 


In  this  case  we  see  that  the  stability  boundaries  are  especially  simple, 
and  only  depend  upon  the  first  two  of  the  six  parameters  UstBd  above. 


The  above  conclusions  suggest  the  following  as  a  reasonable  approach 
to  the  construction  of  stability  boundaries* 

“o 

a.  Plot - —  versus  the  six  parameters  listed  above.  This  mode 

(l)  , 

of  c 

of  presentation  has  the  advantage  that  the  frequency  coalescence 


solution  appears  as  the  straight  line 


I j  thus,  when  the 


of  c 


exact  solutions  deviate  significantly  from  this  line,  we  know 
that  we  are  working  In  a  regime  where  damping  effects  are 


significant.  Since  — 


u 


ofc 


of  c 


,  It  follows  that  the 


“o  Uo 

value  of  —  at  neutral  stability  can  be  determined  from  - 

e  ofc 


provided  that 


ofc 


Is  plotted  separately.  The  latter  plot 


Is  extremely  simple,  since 


ofc 


depends  only  upon  two  param¬ 


eters,  R  and 


n# 


ac 


.  Such  a  plot  Is  presented  In  Figure  26. 


ac 
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Note  that  the  velocity  V  corresponding  to  neutral  stability 
can  be  found  If  the  frequency  of  the  elastic  mode  toe  Is  known, 
since 


a)  »  —  a  and  ca 
owe  o 
e 


2m  r 


Note  that  we  is  not  a  function  of  V,  but  Is  the  frequency  of 
the  elastic  mode  In  vacuo.  If  the  frequency  of  Instability  w 
Is  desired.  It  can  be  found  from  the  equation 


b. 


c  • 


The  three  parameters  R  * 


£  9 

re 


XrUi  "  Xac  jr_ 

x  '  x 

ac  ac 


are 


Important,  and  wide  variation  of  these  parameters  must  be 
taken  Into  account.  We  need  only  consider  negative  values 


of  x  ,  since  otherwise  the  system  Is  statically  unstable, 

AC 


The  ratios  of  aerodynamic  coefficients,  (I  -  2x)  and  — 3.  , 

L. 

will  be  allowed  to  take  on  values  characteristic  of  subsonic 
and  supersonic  conditions.  Thus,  for  subsonic  flow  (Figures 


27  through  35) 

I  -  2x  -  .5 
-  .0625 
CL  -  6.28 


9h 


and  for  supersonic  flow  (Figures  3&  through  Ijij.) 


I  -  2x  =  0 


7P-S  «=  .0625 

a 


1.5 


If  then  allow  x  ■  — ^  to  vary  over  reasonable  ranges, 

dC  C 

we  obtain  practical  ranges  for  the  aerodynamic  stability 
parameters 

I  -  2x 


ac 

C 

*  -  ^ 

ac  a 


I 


d.  The  parameter 


H  x 


Is  not  an  Important  one  for  vehicles 


ac 


moving  In  air,  since  the  mass  parameter 
*  ■  TJsc 

Is  apt  to  be  very  large  and  5?  will  not  be  zero  for  statically 

8C 

stable  vehicles.  Thus,  It  suffices  to  choose  a  single  value 
for  |A  (p  =  100).  The  procedure  of  c  above  then  determines  the 


range  of  the  parameter  ^  ^  ■  ,  Figure  k5  shows  that  the  effect 


ac 


of  varying  — by  a  factor  of  100  produces  less  than  J% 
^  *ac 

“o  xni  "  xac 

variation  In  —  ,  except  at  points  near  —  — •  *  0. 

of c  xac 
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Figure  25.  Flexible  Airframe  With  a  Rigid  Lifting  Surface 
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Figure  26.  Frequency  Coalescence  Solution 
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F I gur«  27.  Frequency  Coalvscvnc*  Solution 
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Figure  26.  Frequency  Coalescence  Solution 


Figure  39»  Frequency  Coalescence  Solution 


igure  3®»  Frequency  Coalescence  Solution 


Figure  5i-  Frequency  Coalescence  Solution 


Flgur*  52.  Frequency  C  oa  r>c«  Solution 


Figure  J3.  Frequency  Coalescence  Solution 


Figure  34*-  Frequency  Coalescence  Solution 


Flgura  37«  Fraquaney  Coal ascanca  Solution 
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Figure  3®*  Frequency  Coalescence  Solution 


Figure  39.  Frequency  Coalescence  Solution 


Figure  Ifi.  Frequency  Coalescence  Solution 


Ilf 


Solution 
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Figure  li3«  Frequency  Coal 


H5 


Figure  Uu  Frequency  Coalescence  Solution 


